Unifying All Classical Spin Models in a Lattice Gauge Theory 
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The partition function of all classical spin models, including all discrete standard statistical models 
and all Abelian discrete lattice gauge theories (LGTs), is expressed as a special instance of the 
partition function of the 4D Z2 LGT. This unies all classical spin models with apparently very 
different features in a single complete model. This result is applied to establish a new method to 
compute the mean-eld theory of Abelian discrete LGTs with d > 4, and to show that computing the 
partition function of the 4D Z2 LGT is computationally hard (#P hard). The 4D Z2 LGT is also 
proved to be approximately complete for Abelian continuous models. The proof uses techniques 
from quantum information. 
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1. Introduction. — Gauge theories describe the most 
fundamental interactions in nature, like QED, weak in- 
teractions and QCD. Only gravity has evaded a quan- 
tum version of the gauge principle. Lattice Gauge The- 
ories (LGTs) are cutoff regulations of gauge theories of 
strongly interacting particles [2, Q ■ When formulated on 
a lattice, LGTs become a new type of statistical mechan- 
ical models that are interesting by themselves, regardless 
of their connection with quantum gauge theories in con- 
tinuous space-time. 

While Standard Statistical Models (SSMs) (like the 
classical Ising or Potts models [3|) have global symmetries 
which arc broken by local order parameters leading to dif- 
ferent phases LGTs have been introduced as models 
with local symmetries [f| which can only be broken by 
global order parameters, corresponding to closed string 
observables. LGTs have also been constructed for arbi- 
trary non-Abelian gauge groups in the context of QCD 
to describe strong interactions [6| . Many interesting phe- 
nomena are known to arise in LGTs which are different 
from SSMs. For example, LGTs may exhibit non-trivial 
phase diagrams that do not correspond to any contin- 
uum field theory. Already the simplest instances of LGTs 
— Abelian discrete LGTs — exhibit remarkable features 
such as a rich phase diagram depending on the gauge 
group 7L q and the underlying lattice, or confinement in 
the strong coupling limit. In fact, it has been argued 
that the center of the group SU(q), i.e. the gauge group 
of Abelian discrete LGTs, Z g , plays an important role in 
the confinement problem 

Given the variety of features of these models, a funda- 
mental question arises: is it possible to give some struc- 
ture to the set of all classical models? In this Letter we 
will give a positive answer to this question by showing 
that an LGT with Z2 gauge symmetry in a four dimen- 
sional square lattice, the 4D Z2 LGT, is complete with 
real coupling strengths for all Abelian discrete LGTs and 



all discrete SSMs. Here completeness means that the par- 
tition function of a large set of models (here, all classical 
spin models) can be expressed as a specific instance of 
the partition function of the complete model (thus, the 
notion of completeness can be seen as a form of universal- 
ity). While similar completeness results for the 2D or 3D 
Ising model have been recently found in the context of 
SSMs d, 0, they are either restricted to specific classes 
of models [8| , or are general but require complex coupling 
strengths, thereby lacking a physical interpretation [9j. 

The completeness results we present here are general 
and are entirely based on real coupling strengths, and 
thus are not only mathematical relations, but have phys- 
ical implications. Note that all completeness results re- 
quire to consider inhomogeneous coupling strengths in 
the complete model. The results are general in the sense 
that they hold for all models with an arbitrary interac- 
tion pattern (including k— body interactions, and in any 
dimension d) between arbitrary q— level spins, which in- 
clude all Abelian discrete LGTs and all discrete (Abelian) 
SSMs. The latter result establishes a general, explicit 
relation between models of very different physical origin 
(e.g. models with local as opposed to global symmetries), 
in contrast with the specific and rather involved connec- 
tions between certain SSMs and certain LGTs known be- 
fore Furthermore, we will also show that the 4D 
I12 LGT can efficiently approximate Abelian continuous 
LGTs with polynomial accuracy, as well as continuous 
SSMs, for which the efficiency depends on the scaling of 
the parameters, as will be specified below. To establish 
the proof, we use results from quantum information. 

2. Completeness of the 4D Z 2 LGT. — In order to prove 
the main result of this Letter, i.e. the completeness of 
the 4D Z2 LGT, we first present a quantum formula- 
tion of the partition function of Abelian discrete LGTs 
(Sect. 2.1), which allows us to construct systematic map- 
pings to other models by properly choosing interaction 
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strengths (Sect. 2.2). Then, we present a method to ob- 
tain general n— body interactions using only Ising-type 
interactions (Sect. 2.3). Finally, we shall use the manip- 
ulations of 2.2 to obtain all interactions required in 2.3 
(Sect. 2.4). 

2.1. Quantum formulation of the partition function of 
Abelian discrete LGTs. — We consider a standard defini- 
tion of an Abelian discrete LGT using a Wilson Hamilton 
function in terms of face interactions, with gauge group 
Wiq. That is, we consider a lattice in d dimensions in 
which each face / G F has k sides or edges. Classical 
spins sit at the edges e G E of the lattice, they have 
q levels, s e — 0, 1, . . . , q — 1, and they are subject to a 
k— body interaction in every face: 
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where J\...k is the interaction strength at face / that has 
spins si, . . . , Sk at its boundary, and the cosine depends 
on the sum of the spins modulo q. For q = 2 (and any k 
and d) we shall refer to these interactions as "Ising-type 
interactions" , and to the model as "d Z 2 LGT" (by de- 
fault meaning k = 4). Note that in this case each term in 
the Hamilton function takes the form Ji...fc(— l) Sl+ "' +Sfc , 
thereby only depending on the parity of the k adjacent 
spins. The Hamilton function (fT]) is invariant under the 
local transformation g v — Yi e e adj o-^ej where e adj v 
are edges adjacent to a vertex v, and X e is defined as 
X e : s e —> (s e + l) mod q . The gauge group is generated 
by these transformations: 7L q = (g v ,Vv € V). The parti- 
tion function of this system is defined as 
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where (3 — l/(ksT) is the inverse temperature, and s := 
(si, . . . , s\e\) is the spin configuration. 

Now we present a quantum formulation of l[2[). In the 
same spirit as in we define a quantum state of 

q— level quantum particles, which we can imagine to sit 
at the center of each face, as in Fig. [D(a), 



I^lgt) = ]T (g) |si 



Sk) 
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where s\,...,Sk are the (values of the) spins at the 
boundary of face / and they are summed modulo q. We 
also define a product state |a) = (^)y gF \ a f) with 

!«/)= J2 e' 3Jl - fcCOS [^ (si+ - +Sfc) ]|si + ... + s fe ). (4) 

The basis states in <j3j> and |4]) are the eigenstates of 
the quantum phase shift operator Z\j) = e l2 ^^ q \j), for 
j = 0, 1, . . . , q— 1. The partition function of the Abelian 
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FIG. 1: (color online) (a) The state [i/'lgt) (here, on a 2D 
square lattice) places one quantum particle (blue dots) at ev- 
ery face, which characterizes the interaction of classical spins 
(black dots) on that face, (b) The left face is merged with the 
right one by setting J1234 = 00 (marked with a shaded face). 



discrete LGT @ is then obtained by computing the over- 
lap between IV'lgt) and \a): 



Zlgt = HV'lgt)- 



(5) 



2.2. Merge and deletion rules. — Similar as in we 
now show how to manipulate the state IV'lgt) by means 
of measurements. That is, we show how to choose a 
product state I7) = <8>y e ^|a/} on a subset F of quantum 
particles such that the state of the remaining particles af- 
ter the projection onto (7I, IV'lgt) = ((tI ® ^)IV'lgt}, is 
again of the form © , but now defined on a different lat- 
tice L'. This allows us to relate the partition function of 
models defined on L' to the partition function of the 4D 
Z2 LGT via (JSJ). More precisely, we consider the merge 
rule of faces, which corresponds to setting \otf>) = \0f>) 
for a certain face /'. From (fj| we see that this choice 
corresponds to setting Jp =00 [9]. Moreover, applying 
this projection on state ([3]) we see that it sets the con- 
dition si + . . . + Sfc = for the k spins at the boundary 
of /'. Now we can substitute this condition on a neigh- 
boring face which, instead of depending on, say, Sk, will 
now depend on s\ + . . . + Sk-i- For example, if we apply 
(0| on the left face of Fig. [Ha), we impose the condition 
si + s 2 + S3 + S4 = 0. This can be substituted in the right 
face, which will now depend on (si+S2+S3)+ss+S6+S7, 
thereby effectively merging the two faces into one larger 
face with an Ising-type interaction with strength J4567 
(Fig. HKb)). We note that the merge rule can be con- 
catenated, thereby enlarging faces at will. On the other 
hand, the deletion rule of faces corresponds to projecting 
a face /' onto the symmetric state J2 Sl Sk ( s i + • • ■ + Sfc | , 
i.e. to setting Jf = 0. 

2.3. Method to obtain arbitrary n—body interactions. — 
Next we show that we can generate a totally general 
interaction between n 2— level particles if all Ising-type 
k— body interactions between these n particles are avail- 
able, for any subset of k particles and all k = 0, 1, . . . , n. 

A general interaction between n spins corresponds to 
assigning a different energy A s to each spin configura- 
tion s. Hence, we need to show that there always ex- 
ists a combination of Ising-type interactions on the dif- 
ferent subsets, i.e. some J,J\,..., J n , J12, . . . , Ji... n (the 
subindex indicates which particles participate in that in- 
teraction) such that J{- 1 ) + Ji ( - 1 ) 81 + . . . + J„ ( - 1 ) s " + 
Ji 2 (-l) Sl+S2 + ... + Ji2...»(-l) Sl+ - +s " = A s is satis- 
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FIG. 2: (color online) Spins fixed by the gauge are marked 
in red. A single-body, a two-body and a three-body Ising- 
type interaction with coupling strengths Ji, J12 and J123 are 
shown in (a), (b) and (c), respectively. 




fied for arbitrary A s and for all s. This is equivalent 
to showing that a matrix with one prefactor of the J's 
in each column (i.e. (— 1)°, (— l) Sl , . . . , (— l) Sn , . . .) and 
with one row for each spin configuration is always invert- 
ible. Notice that this is a square matrix, since there are 
J2l =Q ( n k ) = 2™ J's (columns), and 2™ A s 's (rows). But 
by construction all rows are linearly independent, hence 
it has nonzero determinant, and thus it is invertible. 

2-4- Explicit construction. — We now use the tools of 
2.2 in order to show that we can obtain all interactions 
required in 2.3 in a 4D Z2 LGT. The proof will require 
to fix some spins using the gauge symmetry of the model 
(i.e. fixing them to zero while leaving the Hamilton func- 
tion invariant), a technique which can be applied as long 
as the edges whose spins are fixed form at most a maxi- 
mal tree (i.e. do not form a closed loop) [Til ] . 

First, a "single-body interaction" of si (analogous to a 
magnetic field) is obtained by letting s% interact with all 
other spins around a face fixed by the gauge (Fig. fjja)). 

A two-body interaction is obtained by merging the 
front, lower and back face and creating the face with 
blue boundaries of Fig. d^b). By fixing with the gauge 
six of the spins at the boundary of the blue face, this face 
depends only on si + S2+r + r = s% + S2 (since the sum is 
mod 2) (see Fig.Efb)). Thus, this effectively corresponds 
to a two-body Ising-type interaction between si and S2 
with an interaction strength J12. Notice that by setting 
J 12 = 00 as well, we force si + S2 = 0, which can be seen 
as a propagation of si into S2 (since s% = $2). A concate- 
nated application of this two-body interaction results in 
an effective propagation of a spin through a certain path 
(the turnings of the path can be done similarly). 

A three-body interaction is obtained by bringing three 
spins Si, s 2 , S3 close to each other and then merging the 
large blue face as indicated in Fig.^c). The interaction 
in the blue face corresponds to a three-body Ising-type 
interaction between Si,S2 and S3 (since r% and r 2 are 
summed twice) and with an interaction strength J123. 

The generalization to fc— body interactions with k > 4 
can be done in a similar way as in Fig. ]2jc). The spins 
Sj taking place in the final fc— body interaction are never 
adjacent, and each of them is part of a face at the front, 
back or side with three spins fixed by the gauge (red u- 
shapes) . All but one of the remaining faces are merged, 



FIG. 3: (color online) Replication of spins in four dimensions. 
Yellow faces represent the fourth dimension, and they have 
the same meaning as blue faces, that is, S2 propagates into S3 
by the same method as it propagates into S4. 



and the interaction strength J\,..k in that face determines 
the fc— body Ising-type interaction. 

Thus we have shown how to obtain fc— body Ising- 
type interactions between any group of fc particles, for 
fc = 1,. ,.,n (the zero-body interaction required in 2.3 
is a constant factor, so we obtain Z up to this factor). 
Since the total number of interactions is 2 n , we only need 
to show that a given spin can participate in 2 n /n interac- 
tions. This means that each spin must have this number 
of "end faces" , i.e. faces at the end of a propagation that 
participate in a (many-body) interaction. For example, 
if we use Fig.[^b) to propagate si (i.e. we set J i2 = 00), 
then Si has two end faces, the left and the right one, each 
of which can participate in, say, a three-body interaction 
like the one shown in Fig. [He). But, as can be seen from 
Fig. ]U(b), the propagation of a particle (in 3D) essen- 
tially behaves as a "pipe" which has only two end faces. 
In fact, the number of ends that an encoded particle of 
dimension d e in a lattice of dimension d can have are 
2(d — d e ). Here we essentially have d e = 2, and thus for 
d = 3 the particle is blocked to have only 2 ends. We need 
to resort to a 4D lattice in order to obtain 2(d — d e ) > 2 
ends, and then this replication in different directions can 
be multiply applied until the particle has 2 n /n ends (see 
Fig. [3] for a replication of one spin si into three other 
"end faces" S3, S5 and sq). We refer the reader to [l2| 
for the detailed construction. We remark that all faces 
which are not mentioned in this construction have to be 
deleted using the deletion rule. We also mention that we 
have tried several other procedures to obtain this result 
in 3D, but none of them could avoid the formation of 
loops of edges fixed by the gauge. 

This proves that we can generate a totally general 
n— body interaction between n particles in a 4D Z2 LGT. 
This includes all classical spin models with q = 2 in ar- 
bitrary dimensions d, arbitrary graphs, and arbitrary in- 
teraction pattern. Moreover, by encoding a q— level par- 
ticle in m q — [log q~\ 2— level particles, this also includes 
general interactions between 71' q— level particles, with 
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n' = n/m q . This proves that the 4D Z2 LGT is complete 
for all Abelian discrete classical spin models, including 
all Abelian discrete LGTs and discrete SSMs. 

2.5. Approximate completeness for Abelian continuous 
LGTs and continuous SSMs. — We can go further and 
show that the 4D Z2 LGT is also approximately com- 
plete for Abelian continuous models, that is, the partition 
function of a continuous model can be expressed, up to 
a certain accuracy, as a specific instance of the partition 
function of the 4D Z2 LGT. To see this, we just need 
to let q — > 00 (the lattice spacing remaining discrete) 
and determine what approximation can be obtained (see 
below) . 

2.6. Efficiency results. — The construction presented 
above enables one to generate, from a 4D Z2 LGT, Hamil- 
ton functions that contain M terms with at most k— body 
interactions with an overhead that scales poly(M, 2 k ) for 
q = 2. In the case of q— state models with M general 
A;'— body interaction terms, at most 2 k mq Ising-type in- 
teractions between k'm q 2— level particles are required for 
each term. Therefore, the overhead in the system size of 
the 4D Z2 LGT w.r.t the final model is polynomial if 
k' scales not faster than logarithmically, and q and M 
scale polynomially with the system size. These criteria 
determine whether a given continuous SSM can be ap- 
proximated efficiently. Abelian continuous LGTs usually 
have k fixed (e.g. k = 4), and thus they can be efficiently 
approximated by letting q — > 00 polynomially. 

3. Implications of the main result. — We shall now draw 
three conclusions from the main result. 

3.1. Completeness of the 3D Z 2 LGT with fixed bound- 
ary conditions. — As mentioned above, the only obstacle 
in proving that the 3D Z2 LGT is complete was that 
the replication of spins combined with the k— body in- 
teractions caused loops of spins fixed by the gauge. This 
problem can be overcome by allowing for fixed boundary 
conditions, i.e. fixing spins to zero at the boundary of the 
3D lattice. This means that the completeness results for 
the 4D Z 2 LGT also hold for the 3D Z 2 LGT with fixed 
boundary conditions, which has the same bulk interac- 
tions in the thermodynamic limit as the 3D Z 2 LGT. 

3.2. Mean-field theory for Z 2 LGTs at fixed d > 4. — 
While the mean-field theories of SSMs are easy to con- 
struct, this is not the case for LGTs since Elitzur's theo- 



rem [13j prevents the non-vanishing mean value of a link 
variable. This problem was circumvented using a saddle- 
point approximation 14j, lla ] with the inverse dimension 
1/d as an expansion parameter. The restoration of the 
gauge symmetry is nontrivial in this expansion. Here 
we propose a new method that does not break gauge in- 
variance and that works for fixed d, with d > 4, which 
is based on the construction of a k— clique, i.e. a graph 
with all possible (here, Ising-type) k— body interactions. 
Thus, constructing such a graph is a way of averaging 
over the interaction of a given particle with all the rest, 
and, hence, a way of computing its mean-field theory. 



The construction of the 4-clique for the 4D Z2 LGT has 
been shown in Sec. 2.4 (where we constructed the k— 
cliques for all k = 1, ... ,71). The same construction ap- 
plies trivially to d Z2 LGTs with d > 4, simply by not 
using the extra dimensions (see [13] for further details). 
We mention that this method has potential applications 
in computer simulations. 

3.3. Computational complexity of the 4-D Z2 LGT.— 
Our results imply, in particular, that the partition func- 
tion of the 2D Ising model with magnetic fields can be 
expressed as a specific instance of the partition function 
of the 4D Z2 LGT. Because the former is known to be #P 
complete problem 1|| 17 1 , we conclude that the problem 
of evaluating the partition function of the 4D Z2 LGT 
in the real parameter regime is #P hard, i.e. compu- 
tationally difficult. On the other hand, one can show 
that approximating the partition function of the 3D and 
4D Z2 LGT in a certain complex parameter regime with 
polynomial accuracy is as hard as simulating arbitrary 
quantum computations, i.e. BQP complete [12J. 

4- Outlook. — Our results provide a unification of all 
classical models with very different features into a single 
complete model, the 4D Z2 LGT. In particular, mod- 
els with different types of order parameters, as well as 
models belonging to different universality classes can be 
obtained. It will be interesting to use our results to gain 
further insight in the structure of classical spin models. 
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